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ABSTRACT 

The exo-solar planet HD 806066 has a highly eccentric (e = 0.93) and tight (a = 0.47 AU) orbit. We 
study how it might arrive at such an orbit and how it has avoided being tidally circularized until now. 
The presence of a stellar companion to the host star suggests the possibility that the Kozai mechanism 
and tidal dissipation combined to draw the planet inward well after it formed: Kozai oscillations produce 
periods of extreme eccentricity in the planet orbit, and the tidal dissipation that occurs during these 
periods of small pericentre distances leads to gradual orbital decay. We call this migration mechanism 
the 'Kozai migration'. It requires that the initial planet orbit is highly inclined relative to the binary 
orbit. For a companion at 1000 AU and an initial planet orbit at 5 AU, the minimum relative inclination 
required is ~ 85 deg. We discuss the efficiency of tidal dissipation inferred from the observations of 
exo-planets. Moreover, we investigate possible explanations for the velocity residual (after the motion 
induced by the planet is removed) observed on the host star: a second planet in the system is excluded 
over a large extent of semi-major axis space if Kozai migration is to work, and the tide raised on the 
star by HD 80606b is likely too small in amplitude. Lastly we discuss the relevance of Kozai migration 
for other planetary systems. 



1. INTRODUCTION 

Ongoing radial-velocity surveys have uncovered ~ 100 
exo-solar giant planets. 3 Their orbital characteristics are 
puzzling: many have eccentricities much higher than the 
planetary orbits in our solar system, and a large fraction 
reside very close to the host star. These orbits indicate 
either a formation scenario different from that which op- 
erated in our solar system, or a migration mechanism to 
bring the planets in from large distances. 

The object HD 80606b (Naef, Latham & Mayor et 
al. , 2001, hereafter NLM, period 111.8 days, m p smi = 
3.4Mj) 4 the focus of this article, has the highest eccen- 
tricity (e = 0.927 ±0.012) and the smallest pericentre dis- 
tance (~ 7Rq) among all known exo-planet candidates. Is 
it possible to form such a system in-situ? Among the spec- 
troscopic binaries contained in the Batten catalog (Batten, 
Fletcher & MacCarthy 1989), the only G-type binary that 
has a comparably small periastron distance (~ 6R@) is HD 
27935 (Griffin, Griffin & Gunn et al. , 1985) with e = 0.85. 
Duquennoy, Mayor & Andersen et al. (1992) reported a 
more eccentric system, the K-dwarf binary Gl 586A, with 
e = 0.975 and a periastron distance of IORq. However, 
doubts exist as to whether the latter binary can be pri- 
mordial. But assuming both systems are primordial, their 
presence empirically suggests that HD 80606b might have 
formed as a binary companion to HD 80606. This hypoth- 
esis is currently untestable, the more so since the projected 
mass-ratio in this system is rather extreme. Instead we ask 
the following question: can HD 80606b be formed in the 
cool outer part of a protoplanetary disk and then evolve 
into the current orbit via migration? 

Planet-planet scattering is one such possibility, assum- 
ing HD 80606 had or has multiple planets. Ford, Hav- 



likova & Rasio (2001) integrated the evolution of ~ 10 3 
two-planet systems that are initially dynamically unsta- 
ble. Not a single close encounter produces a planet with 
e > 0.90. This rarity leads us to conclude that either the 
orbits of most observed exo-planets are results of planet- 
planet scattering, or that the HD 80606 system is not the 
result of such scattering. 

The leading hypothesis for planet migration involves 
gravitational interactions between the planets and the gas 
disks out of which they form (Goldreich & Tremaine 1980; 
Lin, Bodenhcimer & Richardson 1996). Under certain cir- 
cumstances planet-disk interactions may also excite the or- 
bital eccentricity of the planet (see, e.g., Goldreich & Sari 
2002). If so, this may explain the orbital separation and 
eccentricity of the bulk of the systems discovered so far. 
These systems exhibit an eccentricity distribution that is 
roughly flat below e ~ 0.6 and distinctly drops off around 
0.6. It is believed that Jupiter-mass planets can open up 
gaps in their natal gas disks with a fractional width Aa/a 
reaching up to a similar value. This coincidence suggests 
that passage of an eccentric planet through the disk on 
either side of the gap tends to damp the planet's orbital 
eccentricity and to limit the maximum eccentricity planet- 
disk interaction produces to e ~ 0.6. However, HD 806066 
stands out with e = 0.93. 

The origin of the orbit of HD 80606b becomes more 
perplexing if we consider tidal effects. Assuming that the 
planet has a tidal dissipation efficiency similar to other 
known exo-planets (quality factor Q p — 3 x 10 5 , Wu 2002) 
or to Jupiter (Goldreich & Soter 1966, Peale & Green- 
berg 1980) and setting its true mass equal to the minimum 
mass, its radius to a Jupiter radius, we find that a mere 
650 Myrs ago, HD 80606b had an implausibly high eccen- 
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4 The host star HD 80606 is a G5 star with [Fe/H] = 0.43, m B = 9.65 and m v = 8.93. 
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tricity of e = 0.99. 5 The age of the system is likely to be 
this old or older; the stellar projected rotational velocity 
is low, vsini w 0.9km/ s, typical of old G-dwarves; it is 
also chromospherically quiet. 

Interestingly, there is a neighbor to this system, HD 
80607, a main-sequence companion ~ 1000 AU away. This 
prompts us to develop a theory in which the companion is 
responsible for the high eccentricity and small orbit of the 
planet. There are two ingredients in our theory. First, we 
assume that the planet was born on an orbit of a few AU, 
and had an orbital plane inclined relative to the stellar bi- 
nary plane. The remote stellar companion would induce an 
eccentricity oscillation in the planet's orbit via the Kozai 
mechanism (Kozail962; see also discussions by Holman, 
Touma & Tremaine 1997, Innanen, Zheng, & Mikkola et 
al. , 1997, and Mazch, Krymolowski & Rosenfeld 1997). 
The second component is tidal circularization. The tides 
operate most effectively during episodes of high eccentric- 
ity in the Kozai cycle. We rely on the dissipative tidal 
process to irreversibly draw the planet inward. Combining 
the two processes, we can explain the current high eccen- 
tricity as well as the tight orbit. We refer to this planet 
migration scenario as 'Kozai migration'. 

Eggleton, Kiseleva & Hut (1998), Kiseleva, Eggleton & 
Mikola (1998), and Eggleton & Kiseleva-Eggleton (2001, 
EKE hereafter), have developed this scenario to explain 
some hierarchical triple star systems. Blaes, Lee & 
Socrates (2002) applied an analogous scheme (gravita- 
tional radiation being the dissipation mechanism) to hi- 
erarchical triple black holes in galactic centers. In this ar- 
ticle, we adopt the formalism of EKE in an attempt to pro- 
duce a plausible life-history for HD 80606b (§2). The story 
depends on the value of the tidal quality factor Q (§3.1). 
The observational consequences of such a life-history fol- 
low in §3.2. In §3.3 we discuss the tidal velocities induced 
at the stellar photosphere by the highly eccentric planet, 
and the possibility of observing such tides. NLM note the 
presence of substantial residuals to the Keplerian fit in 
the HD 80606 system. The tidal velocities appear to be 
large, but not large enough to explain the observed veloc- 
ity residuals. Lastly (§3.4), we assess the importance of 
Kozai migration for other planetary systems where binary 
companions are also known to exist. 

2. KOZAI MIGRATION FOR HD 80606B 

The host star of HD 80606b is known to be a member of 
a common proper-motion binary with a companion (HD 
80607) that is similar to HD 80606 in both spectral type 
and mctallicity . The two stars are separated by ~ 20" on 
the sky. Unfortunately, it is hard to translate this into a 
linear dimension. The Hipparcos distances for both stars 
are highly uncertain 6 This complicates an age determina- centricity. The planetary spin is quickly (pseudo- 
tion with isochrone fitting. But if we enforce the constraint Synchronized with the orbital motion, 10 while the star's 
that the stars have reached the main-sequence, isochrone synchronization time greatly exceeds the orbital circular- 

5 This problem is further exasperated if the planet has acquired the same pcriastron distance while it was still hot and large. 
6 Hipparcos parallaxes are 17±6 mas for HD 80606 and 10±9 for HD 80707. The large error bars are presumably due to mutual contamination 
(see also NLM). 

7 We denote quantities related to the planet and the companion star with subscripts p and c, respectively. 

s Below this value, the pcriapsc argument of the planet circulates through all angles very fast and the torque from the companion is effectively 
averaged to zero. There is little variation in the planet's inclination and eccentricity. Also see, e.g., Holman et al. (1997). 
9 See Ford, Kozinsky & Rasio 2000 for a discussion of higher order effects). 

10 Pesudo-synchronization applies to an eccentric system where the planet spin is tidally synchronized to a rate that is in-between the orbital 
frequency and the angular frequency at periapse. See, e.g., Hut (1981). 



fitting yields a lower limit of ~ 65 pc on the distance, im- 
plying a projected binary separation of 1300 AU at the 
current epoch. On the other extreme, requiring that the 
stars have not yet ascended the red giant branch yields 
a maximum distance of ~ 125 pc, and a projected sep- 
aration of 2500 AU. As noted above, the low vsini and 
the chromosphcric quietness of HD 80606 are typical of an 
old G-dwarf. In this article, we adopt a value of 1000 AU 
for the binary semi-major axis (o c ). The eccentricity is 
taken to be e c — 0.5. 7 We discuss the effects on the Kozai 
migration when a c is increased. 

Even at such a large distance, the companion star could 
significantly perturb the planet orbit as long as the two or- 
bital planes are initially inclined to each other more than 
39.2 dcg. 8 Secular effects of the Kozai type could then 
occur and produce large cyclic variations in the planet's 
eccentricity (e p ) and relative inclination, as a result of an- 
gular momentum exchange with the companion orbit. Ef- 
fects on the companion orbit produced by the planet can 
be neglected, given the extreme ratio in the two orbits' an- 
gular momenta (a c ~> a p and m c 3> m p ). Under this ap- 
proximation, the z-component of the planet's angular mo- 
mentum, J z = GM±apM p ^Jl — e% cos I, is conserved. 9 

Here the z-axis is normal to the companion plane. As a p 
is not modified by secular perturbations, the Kozai integral 

Hk = ^/l — ep cos/ is conserved during the oscillations. 

Minima in e p concur with maxima in /, and vice versa. 
Each Kozai cycle lasts a time ~ P c 2 /P p , or a few Myrs. 

A much slower process, tidal circularization, gradually 
removes energy from the orbit and draws the planet in- 
ward. An initially more highly inclined orbit (larger I) 
could reach higher eccentricity. And since the dissipation 
process depends sensitively on the nearest approach dis- 
tance between HD 80606b and its host star, such an orbit 
would suffer stronger orbital decay and have shorter mi- 
gration timescale. The key issue in Kozai migration is 
therefore the initial value of I. In §2.1, we first estimate 
the minimum initial inclination required for migrating HD 
80606b from an initial orbit of a p — 5 AU and an ini- 
tial eccentricity of e p — 0.1, we then confirm this estimate 
numerically in §2.2. 



2.1. Minimum Initial Inclination: heuristic argument 

To estimate the minimum initial inclination, we note 
that while the Kozai oscillation roughly conserves yfil — 
Cp)cosI and a p during individual cycles, tidal circu- 
larization preserves orbital angular momentum (J = 

V GM± M p «/b p (1 — ep)) during episodes of maximum ec- 
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Symbol 



Definition 



Fiducial Values 



stellar (host & companion) and planetary masses 
stellar and planetary radii 
tidal Love number 
gyro-radius, moment of inertia = r g MR 2 
tidal dissipation quality factors 
spin frequency 

planet semi-major axis, eccentricity, mean motion & period 
parameters for the companion orbit 



M*, M c , M p 



^* i Qp 

dpi ^"pi Tip Pp 

a c , e c , n c & P c 



1.1M , 1.1M , 7.80Mj 
i? Q , Rj 
0.028, 0.51 
0.08, 0.25 
10 6 , 3 x 10 5 
initially 20 days & 10 hours 

n p = ^G(M* + M p )/a p 

a c = 1000 AU, e c = 0.5 



ization time. Hence the planet's orbital angular momen- 
tum cannot be absorbed into spin. Since e p ~ 1 near 
maximum, conservation of J translates into a constant pe- 
riastron distance between maximum e p during the Kozai 
cycles. 

In an orbit with a semi-major axis of 5 AU, HD 80606b 
would have to reach an eccentricity as high as 0.993 (at 
some point during the Kozai cycle) in order to produce 
the currently observed periastron. Adopting a minimum 
of / m 40dcg during the Kozai cycles (see, e.g., Holman 
et al. , 1997), the Kozai integral ( \f(l — e 2 ) cos/) remain- 
ing constant yields I > 84.8 deg for an initial e p = 0.1. In 
other words, the two orbits had to be virtually perpendicu- 
lar to each other when the planet formed. This conclusion 
is independent of the companion orbital separation. 

The planet in HD 80606b is not currently undergoing 
Kozai oscillations; the Kozai mechanism operated only 
when the planet had a semi-major axis a > 4 AU. To see 
why note that the torque exerted by the companion star 
on the planet has its value and direction dependent on the 
pericentre argument of the planet (w p ). If, however, this 
argument also precesses under other forces, the averaged 
Kozai torque is reduced. General relativistic effects, tidal 
effects, and rotational quadrapolar bulges on the planet 11 
are mainly responsible for these extra precession. Their 
rates are summarized here (Sterne 1939, Einstein 1916), 
with definitions for the symbols listed in Table 1., 
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Taking a Kozai precession rate of ~ 0.5n 2 Jn p /{l - elf 2 
(Holman et al. , 1997), and assuming that the planet is 
pseudo-synchronously spinning with the orbit (e p = 0.1), 
the relative precession rates are, 
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If any of the above ratios exceeds unity, the Kozai oscilla- 
tion is destroyed. Setting the first ratio to unity, we obtain 
a suppression radius w 4 AU. Currently, the Kozai oscil- 
lations are strongly suppressed by these extra precessions. 
This conclusion does not affect the minimum inclination 
required for Kozai migration. However, the suppression 
radius does depend on the companion separation: e.g., it 
is reduced to s=s 1.5AU if the companion semi-major axis 
is as small as 250 AU. 



2.2. Numerical Integration 

We adopt the set of equations (eq. [11]-[17]) in EKE 
(also see Eggleton, Kiseleva & Hut 1998) to describe the 
secular evolution of the planet orbit, as well as the spin of 
the host star and the planet. Any change in the compan- 
ion's orbit is ignored. Physical effects described by these 
equations include: secular interactions, GR, tidal and ro- 
tational precession, and tidal dissipation. These equations 
are explicitly listed in Appendix A. 

We integrate equations (A1)-(A7) starting from the 
following initial conditions: planet orbit a p = 5.0 AU, 
e p = 0.1, relative inclination / = 85.6 deg (slightly higher 
than the 84.8 deg estimate in §2.1 because the actual mini- 
mum inclination is ~ 48 deg, not <~ 40 deg as we assumed) , 
and periapsc angle lu p — 45 deg. 

Values of the system parameters are listed in Table 1. 
The resulting orbital and spin information are presented 
in Fig. 1 as functions of time. Initially, the planet under- 
goes Kozai oscillation with each cycle lasting ~ 20 Myrs. 
Whenever e p evolves close to maximum, a small amount of 
energy (but not angular momentum) is removed from the 
orbit. Both a p and e p are decreased, while a p (l — e p ) re- 
mains roughly constant. This corresponds to a small kick 
to reduce e p whenever lu p w tt/2 (see Fig. 3 of Holman 
et al. 1997). The Kozai integral slowly increases. This ex- 



plains why the amplitude of e p oscillation shrink as the 
planet migrates inward (top-left panel of Fig. 1). Af- 
ter about 0.7 Gyrs, the planet has reached an orbit with 
a p ~ 2.5 AU, and the Kozai oscillation is completely de- 
stroyed by GR precession. At this epoch, we find e p ~ 0.99 
and / ~ 50 deg. From this point on, tidal circularization 
dominates the evolution, and the influence of the binary 
companion is negligible. The planet reaches an orbit of 
e p = 0.927 and a p = 0.47 AU - 1.2 Gyrs after its birth. 
A few notes are in order. 

We adopt a value of Q p = 3 x 10 5 for the planet, and 
Qi, = 10 6 for the host star. We discuss the reasons for 
these choices in §3.1. Using them, the planet's spin reaches 
synchronicity with the orbital frequency near periastron 
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almost instantly, while the star's spin is hardly affected by 
tides (see lower-right plot of Fig. 1). Tidal circulariza- 
tion is dominated by dissipation in the planet over that 
in the star by a factor of ~ 10. Most of the dissipation, 
which conserves orbital angular momentum, occurs when 
the planet has been driven to maximum eccentricity by 
the companion star. 

The results do not depend on the value for the longitude 
of the ascending node. This is because in the quadrapole 
approximation for secular interactions, the companion or- 
bit is effectively taken to be circular so there is no preferred 
axis in the binary plane. This also explains why the verti- 
cal angular momentum of the planet orbit (Kozai integral) 
is conserved - there is a rotational invariance with respect 
to the z axis. 

The results also do not depend on the value for the argu- 
ment of pericentre, uo p . Kozai oscillations can cause either 
circulation or libration in the value of u) p . Fig. 3 of Hol- 
man et al. (1997) shows that the circulating solutions can 
reach higher e p and therefore should lead to shorter evo- 
lutionary timescale. However, for the initial conditions we 
adopted, the circulating and librating solutions essentially 
coincide and give rise to similar values for the maximum 
eccentricity. 

To migrate the planet from a distance larger than 5AU, 
one requires a higher minimum inclination so as to reach 
the same pericentre distance at maximum eccentricity. In 
comparison, this inclination hardly varies when the initial 
eccentricity of the planet is varied from 0.1 to 0.6, this 
can be seen from Fig. 1. Lastly, if the companion is fur- 
ther away than 1000 AU, each Kozai cycle lasts longer, 
and it takes a longer time for the planet to migrate to its 
current location. For instance, placing the companion at 
a c = 1500 AU, we find the planet should now be ~ 5Gyrs 
old. 

3. DISCUSSION 



3.1. Q Value 

A time-dependent tide raised on a celestial body, due to 
cither its asynchronous rotation (asynchronous tide) or its 
eccentric movement (eccentricity tide), can be dissipated, 
leading to synchronous spin or orbital circularization. The 
rate of dissipation is conventionally described by a dimen- 
sionless quality factor Q, the ratio between the energy in 
the tidal bulge and the energy dissipated per orbital pe- 
riod. For equilibrium tides, Q is related to the lag-angle 
between the tide-inducing object and the tidal bulge 5 by 
S w 1/2Q. 

It is possible to adapt the Q description to dynamical 
tides, in which case the tidal energy refers to the energy 
in the (non-existent) equilibrium tidal bulge, and Q is ob- 
tained from the dissipation averaged over a range of tidal 
frequencies. Back-reaction to the orbits by the gravita- 
tional moments of the tidal waves is ignored, a proce- 
dure strictly justifiable only when the damping time of the 
waves is shorter than the orbital period. However, secular 
effects of the tidal dissipation should be independent of 
this back-reaction. 

12 This statement is independent of the Q* value of the host stars. The 
effects of both circularizing and eroding the planet orbit. This would fail 
dissipation in the planets. 

13 If this possibility turns out true, tidal dissipation in stars will be held accountable for the observed inner cutoff in the orbits of exo-planets. 



Jupiter's Q value has been estimated to be 10 5 <J Q <J 
2 x 10 6 , with the actual value believed to be closer to the 
lower limit (Goldreich & Soter 1966, Peale & Greenberg 
1980), based on the resonant configuration of the Galilean 
satellites. The physical origin of this Q value has remained 
elusive for a few decades. As is explicitly shown in Wu 
(2002), exo-planets share similar Q values (Q p ~3x 10 5 ) 
as Jupiter. 12 This is striking as the exo-solar planets have 
different thermal environments, formation histories, and 
possibly interior compositions, than those of Jupiter. The 
similar Q values provide an important clue as to the under- 
lying dissipation mechanism. Moreover, while the Q factor 
for Jupiter pertains to the asynchronous tide, the one for 
exo-planets concerns the eccentricity tide. One may spec- 
ulate that the quality factors for the different Jovian-mass 
objects agree because all the objects rotate fast with peri- 
ods comparable to or shorter than the tidal period. 

Tidal dissipation in the host star of HD 80606b also 
contributes to circularization. What is the tidal Q factor 
appropriate for solar-type stars? Field solar-type binaries 
are observed to be circularized out to <~ 12 day orbits 
(Duquennoy & Mayor 1991). This would imply Q ~ 10 5 
adopting an age of 5 Gyrs. These binary stars have pre- 
sumably been synchronized, so this Q value is affiliated 
with dissipation of the eccentricity tide in rapidly spin- 
ning objects (Qfast)- For the evolution of planet orbits, 
the Q value of relevance is that for the asynchronous tide 
in slowly rotating objects to which most host stars belong 
(Qslow)- It is not clear that Qf ast is directly related to 
Qslow What observational constraints could we put on 

Short period planets can tidally spin up the star, at the 
cost of planetary orbital decay. Assuming all host stars 
are Sun-like and all exo-planets are Jupiter-like, we obtain 
a lower limit of Q s iow > 10 5 if all planets are to survive in 
their current orbit for another 5 Gyrs. 

Drake et al. (1998) considered the rotational velocities 
for three planet-bearing stars (r-Boo, u-And, and 51 Peg). 
Among them, r-Boo has the closest and most massive 
planet and appears to be synchronized, while the other 
two have not. Taken at face value, this would imply 
Qsiow ~ 10 4 . This contradicts our earlier statement that 
Qsiow > 10 5 . These two results can be reconciled either 
by the imminent death of the shortest period planets, 13 or 
by the possibility that r-Boo has had only its surface con- 
vective layer synchronized with the orbit (see Goldreich & 
Nicholson 1989 for a similar argument on massive stars). 
The latter possibility also makes it difficult to infer Q s i ow 
using the spin data of solar-type binaries. It is interesting 
to notice that among the three stars discussed by Drake 
et al. (1998), r-Boo likely has the thinnest convective en- 
velope. 

Unable to infer hard constraints on Q s iow, we decide to 
adopt Q± = Q s iow = 10 6 in our study. With this choice, 
dissipation in the planet is roughly 10 times more impor- 
tant than dissipation in the star even for a planet as mas- 
sive as M p = 7.8Mj. This choice has the advantage that 
unless the actual is much lower, the results presented 
here are not qualitatively affected. 

stars are typically slowly rotating so their tidal dissipation have the 
to produce the observed circular orbits, if it is more important than 
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Fig. 1. — Evolution of the planetary orbit as a function of logarithmic time (in years). The four panels represent (from left to right, up to 
down), evolution of the planet eccentricity, planet semi-major axis and periastron distance (in AU), inclination of the planet plane relative 
to the stellar binary plane, and spin periods for both the planet and the star (in days). See the text for the choice of parameters. During 
the Kozai cycles, maxima in eccentricity concur with minima in inclination. At the current orbit (e p = 0.927 and a p = 0.47 AU), the planet 



no longer undergoes Kozai oscillations as a result of GR and tidal precession, 
planet orbit under tidal circularization. 



The dashed part of each curve depicts future evolution of the 



3.2. Constraints on possible Second Planets 

If there existed a second planet in the system, secu- 
lar interaction 14 between this planet and HD 80606b may 
destroy the Kozai oscillations. Kozai oscillations require 
that the precession caused by a second planet acting on 
HD 80606b's orbit be no larger than that produced by the 
companion star (Innanen et al. , 1997). We briefly investi- 
gate what constraints this puts on the mass and orbit of 
the second planet. 

We assume the two planets have coplanar orbits and 
the second planet (denoted with subscript s) lies inward 
of the orbit of HD 80606b. The derivation is similar if it 
lies outward. Secular interaction between the two plan- 
ets alone gives rise to variations in their eccentricities (e) 
and periapse angles (oj). Adopting a complex variable 
I = e exp(zwi), and assuming that e< 1, the first-order 
secular evolution equations read (Murray & Dermott 2000, 
hereafter MD2000), 



dl 



dL, 



^ iA ss I s -\- iA S pIp : — iAp S I s -\- iAppIp : (3) 



where the real coefficients An are 



A lt 



a 2 M, 

a M s 
~4M* 



a 2 M n 



-n s Bi, A sp = — -^r-n s B 2 



n p B 2 , 



Here, n s = y/ G(M* + M s )/af 



$2 («). B 2 



4 M* 
. a M s . . 

A pp = 1Jf n P Bl - ( 4 ) 

a = a s /a p , and B\ = 



^3/2 ( a ) with bp being the usual Laplace 



coefficients (MD2000). The general solution to eq. (3) is 
composed of two linear eigenvectors with their respective 
precession rates 



{A ss + App) ± (A ss — App) 2 + AA S pAp S 
9± = o ■ ( 5 ) 



We are able to approximate the precession rate of HD 
80606b under secular interaction as, 



dujp 
~dt 



= Min( 5+ ,.g_) ifM p >M s , 

Max( 5+ ,g_) ifM p <M s . (6) 



At birth, HD 80606b could not have undergone Kozai 
oscillation if the second planet is massive enough. Requir- 
ing the above secular precession rate to be smaller than 
the Kozai precession rate (see eq. [2]) when a p = 5 AU, 
we find that the mass of the second planet is constrained to 
a value that depends on its distance to HD 80606b. This 
is depicted in Fig. 2. 

These results also strongly constrain the current state of 
the planetary system. No other Jupiter-mass planet could 
now live between about 0.05 and 100 AU. Even earth-mass 
cores are excluded between 1 and 20 AU. Did such objects 
never form? Or were they cleared away? We do not ex- 
plore these issues here. 

In any case, Kozai migration is incompatible with the 
notion that the claimed residual velocity in HD80606 
(NLM) is due to a second planet in the system. 



We ignore interactions of the mean-motion type as they are important only for a small phase space. 
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a s (AU) 

Fig. 2. — The constraints on the mass (vertical axis, in unit of Mj) and orbital semi-major axis (horizontal axis, in AU) of the possible 
second planet placed by the Kozai migration of HD 80606b. The shaded region must be cleared before HD 80606b (assumed to be at 5 
AU) can undergo Kozai oscillations. The latter force HD 80606b into a variable but highly eccentric orbit with the minimum periastron and 
maximum apastron delineated by the two dotted vertical lines (see Fig. 1). No planets should have survived within this region. Planets lying 
above the slanted dashed line would incur a radial velocity larger than 20m/s on the host star (taking sini = 0.5) and should have been 
detected were they present. 



3.3. Can the Tide Raised by the Planet Account for the 
Residual Velocities? 

NLM reported residual velocities of order ~ 40 m/ s 
in the spectra of HD 80606. These are higher than 
the expected noise and there is no clear trend in time. 
Could they be caused by anything other than an unknown 
planet? We consider one option: the tidal velocities in- 
duced by HD 80606b on the surface of the star. The stellar 
rotation is ignored in the following analysis. 

First, we assume that the star adjusts its hydrostatic 
equilibrium instantaneously according to the tidal poten- 
tial of the planet at a radial separation D = D(t) (the 
equilibrium tide picture), and estimate the corresponding 
tidal velocities. The maximum radial displacement at the 
stellar surface is roughly 

with the horizontal displacement being of a similar order. 
Velocities associated with these displacements are simply 
estimated as v cqu i = -l^ 2 -, yielding a maximum linc-of- 
sight velocity ~ 60cm/s occurring within 10 hours of the 
closest-encounter. This maximum requires conjunctions to 
coincide roughly with the line of periapse, as is the case in 
HD 80606. 

We turn next to consider a more realistic description 
of the tidal amplitudes. The stellar response toward the 
time- varying potential is decomposed into a series of eigen- 
mode oscillations, with the mode amplitudes maximized 
when the tidal forcing frequencies are in resonance with 
the mode frequencies (the dynamical tide picture). If the 



planetary orbit is nearly circular, the dominant forcing 
frequency is ~ 2n p . At an orbital period of 111 days, 
this forcing resonantly excites gravity-modes of radial or- 
der ~ 1000. The spatial overlap between the tidal po- 
tential and such a high-order mode is too weak to be 
interesting. However, HD 806066 has a highly eccentric 
orbit. Tidal excitation occurs primarily near periastron 
and the dominant forcing frequencies is about twice the 
orbital frequency near periastron, and is related to Tip 3jS 
<~ 2n p /(l — e p ) 3 / 2 <~ 100rip. Gravity- modes of radial orders 
^ 30 — 80 become relevant and their resonant excitation 
give rise to much larger tidal velocities. A brief treatment 
of the dynamical tides is given in the appendix (§B). Here 
we summarize our results. 

We adopt the standard solar model (Christensen- 
Dalsgaard 1998) to calculate the eigenmode frequencies 
and damping rates. The response of the star depends sen- 
sitively on how close an orbital multiple lies around one 
such eigenmodes. The maximum response occurs when the 
frequency off-resonance (5a) is of order or smaller than the 
mode damping rate (7), giving rise to velocities of order 
10 to 100 m/ s. However, this is a rare event. In the fre- 
quency range of interest (150 to 550n p ), there are roughly 
Npe&k ~ 200 discreet forcing peaks, and -/V mo d ~ 50 
eigenmodes. These modes typically have line-widths cal- 
culated from radiative diffusion to be 7„/cr„ ~ 10~ 7 . A 
typical best resonance 5a = Min(|cr„ — aj\) is of order 
<~ « P /A mo do ~ cr„/10 4 <~ 10 3 7„. So the chance of one 
mode lying close to an orbital multiple with a n < 7„ 
is roughly 10~ 3 . Typical tidal velocities therefore range 
from lem/s to lOcm/s, well below the level of the ob- 
served residues. One such example is shown in Fig. 3. 
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Fig. 3. — Dynamic response of a Sun-like star under tidal forcing from a planet of mass 3.9/ siniMj, an orbit of period 111.82 days and 
eccentricity 0.927. We only consider the dominant I = \m\ = 2 component in the tidal forcing. The upper panel plots mode amplitude (in 
terms of surface horizontal velocity) as a function of mode period in days. Usually one or two modes dominate the response, and in this case 
it is a n = 41 gravity- mode. The dotted and dashed lines are the frequency spectrum of tidal forcing over one period (Kj , see appendix), 
and the overlap integral between the tidal potential and the eigenmode (Q„e), both in arbitrary units. The forcing spectrum peaks around 
2(1 — e p ) — 3 ' 2 n p and the overlap integral decreases sharply when mode period (and mode radial order) increases, roughly as Q n t oc n~ 5 . 
Besides from these two functions, the tidal response also depends on the goodness of the resonance and the mode damping rate, these are 
plotted in the lower panel. The frequency detuning 5a = \a n — Cj\. Here <r n is the frequency of a radial order n mode, and <Tj is the closest 
integer multiple of the orbital mean motion (<Tj = jn p and j is an even number). The damping rate 7„ is estimated from radiative diffusion. 
The two rows of short vertical bars illustrate the positions of the crj (upper row) and a n (lower). 



The highest amplitude mode has a fractional Lagrangian 
pressure perturbation Sp/p ~4x 10~ 7 near the surface. 
Translating this into light variation, this mode will be as 
hard to detect as the solar oscillation. 

To showcase the sensitivity of the tidal responses, we 
plot in Fig. 4 the marked variations in the tidal responses 
when the orbital period is tuned through the observed 
uncertainty range: 111.81 ± 0.23 days (NLM). Slight dif- 
ferences in the stellar structure affect the eigenmode fre- 
quencies and would also give rise to similar variations in 
the tidal responses. So the response results are best ap- 
preciated as a probability distribution of amplitudes, also 
shown in Fig. 4. We find that there is only a <~ 1CP 3 
chance that the tidal velocity reaches the claimed level of 
residues, ~ 40 m/ s. 

We calculate the damping rates taking only radiative 
diffusion into account. Turbulent viscosity may enhance 
this damping. But this will not appreciably change our 
conclusion: raising the damping rates by a factor of 10, 
we obtain similar results as in Fig. 4 except that high 
tidal amplitudes are less likely to occur. 

On the other hand, when we substitute the solar model 
for a 1.25M Q ZAMS model (courtesy of Christensen- 
Dalsgaard), we found that the typical tidal amplitudes 
rise by about a factor of 10. Equivalcntly, high velocity 
amplitudes appear roughly 10 times more frequent. Com- 
pared to a solar model, the 1.25M Q model has a thinner 
convection zone. Its gravity-modes have a relatively larger 



surface amplitude. 

In summary, thanks to the high eccentricity of HD 
80606b, the tidal velocities reach tantalizingly close to the 
current detection thresholds, but likely fail to explain any 
residual velocity of order 10 m/ s. However, if the residual 
velocity is indeed associated with the dynamical tide, we 
expect it to be a coherent oscillation throughout the orbit 
and have a period in the range of 0.2 to 0.6 days. The 
inclination of the orbit (sini) does not affect our conclu- 
sions. 

3.4. Other Planetary Systems 

Zucker & Mazeh (2002) showed that planets in binaries 
are statistically different from those around isolated stars 
in their period-mass distribution. The former population 
dominates the lower-right quadrant of Fig. 5. This dis- 
tinction suggests that planet migration in binaries may be 
induced differently from those around single stars, Kozai 
migration being one such possibility. Among the circled 
planets in Fig. 5, we have only studied r-Boo and 16 
CygB. We found that t-Boo, the most massive close-in 
planet, could have gone through Kozai migration and later 
circularized, while there is no compelling evidence of Kozai 
migration for 16 CygB, though it could currently be un- 
dergoing Kozai oscillations (Holman et al. (1997). A bi- 
nary companion may assist planet migration through other 
means, e.g., exciting tidal waves in the proto-planetary 
disk and affecting angular momentum transport. These 
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Fig. 4. — The probability distribution of the tidal response. The upper panel plots the horizontal velocity on the stellar surface as a function 
of orbital period, with the period running through the observational uncertainty range. Sharp spikes occur when an eigenmode is resonantly 
forced (<5<r < J n ). Other parameters are the same as in Fig. 3. From these data, one finds a probability distribution for the tidal velocity to 
surpass a certain value. This is depicted in the lower panel. Typical amplitude (probability of 0.5) is ~ 20cm/s, with the claimed level of 
residue, ~ 40 m/s, occurring at a probability of ~ 10~ 3 . 



deserve further study. 

Could some of the planets around currently single stars 
have undergone Kozai migration followed by the subse- 
quent disruption of the binary? We consider this unlikely. 
The rate of Kozai migration is limited by the tidal circular- 
ization process. For reasonable values of closest periapse 
distance, the latter process takes a few 10 8 years or longer. 
By this age a typical stellar cluster has dispersed and the 
rate of stellar encounter has become negligible. However, 
it is possible that a planet's eccentricity was excited by 
Kozai interaction with the binary companion before the 
latter was removed from the system. A simple estimate, 
however, finds that this may account for the observed ec- 
centricity in at most <~ 10% of the known systems. 

4. SUMMARY 

It is striking that a companion star as remote as 1000 
AU can force a planet to migrate inward. Kozai migration 
is a long timescale process and demands a 'clean' planetary 
system in which no other planetary object exists. Kozai 
migration can be responsible for the present orbit of HD 
80606b only if the initial relative inclination is sufficiently 
high (I > 85deg). Assuming a population of similar triple 
systems with randomly distributed /, ~ 7% of them could 
have gone through Kozai migration. So on the strength of 
one object (HD 80606b), one may speculate about the high 
density of triple stellar-planetary systems yet undetected. 

Our prediction regarding the lack of a second planet in 
the system (Fig. 2) should be tested. In this regard, we 
note the claim of the abnormally large residue velocities 



on the star. We investigated the possible association be- 
tween these velocities and surface tidal waves excited by 
the highly eccentric planet. It seems unlikely that the lat- 
ter can account for the former if the star has a structure 
like that of our Sun. And if the system is indeed devoid of 
second planet, it raises issues as to how it is cleaned. 

Zucker & Mazeh (2002) pointed out that planets in 
binary systems and planets around single stars do not 
share the same period-mass distribution. The implications 
might be that binaries can cause planet migration - Kozai 
migration being one of the more definite possibilities. 

HD 80606b's highly elongated orbit and strong tidal in- 
teraction combine to make it a possible target for direct 
detection. Its current tidal luminosity (which is likely to 
be converted into heat and radiated all through the or- 
bit) is ~ 10 28 erg/ s, or ~ 14 mag dimmer than the parent 
star, at a maximum angular separation of 0.01" (assuming 
distance of 80 pc). A comparable contribution may arise 
if the planet evenly emits the stellar insolation it receives 
during periastron passages. 

Part of the research was performed while Y.W. was on 
a CITA-ITP postdoctoral exchange programme funded by 
NSF grant PHY 99-07949. Y.W. would like to thank Pe- 
ter Eggleton for his patience, Ludmila Kiseleva-Eggleton, 
Omer Blaes and Man-Hoi Lee for useful conversations. 
N.M.'s research is supported by NSERC of Canada and 
by the Canada Research Chair program. This research 
benefited from ADS and SIMBAD. 
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Fig. 5. — Period versus projected mass for all observed planet candidates (http://cfa-www.harvard.edu/planets/, Sep. 2002). Circles are 
drawn around planets in visual binaries (Worley & Douglass 1997) with a projected separation < 20". Besides from the 9 binaries mentioned 
in Zucker & Mazeh (2002, HD 142, Gl 86, HD 19994, e Eri, r Boo, HD 178911B, 16 Cyg B, HD 195019, HD 217107), there are also HD 80606, 
HD 213240, HD 121504, HD 46375, HD 27442, BD-103166, and 7-Cepheid. Planets lying below the dashed line incur velocity variation higher 
than 20m/s. There appears to be an over-abundance of close-in, massive planets that are also in binaries. 
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APPENDIX 
EQUATIONS FROM EKE 

We adopt the set of equations (cq. [11]-[17]) in EKE (also see Eggleton, Kiseleva & Hut 1998) to describe the secular 
evolution in the planet orbit, the spin of the host star and the planet. The companion is assumed to be unaffected. 
Physical effects described by these equations include: secular interactions, GR, tidal and rotational precessions, and tidal 
dissipation. These equations are explicitly listed here, 

V*-V p - 5(1 - e 2 p )S eq , (Al) 
2e 2 

2W*-2W p -—^(V ir + V p ), (A2) 

1 Cp 



1 de p 
e p dt 

1 da p 
a p dt 
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dl 
~dl 

dcj p r „ . „ „ ,„,„ „ sn cos/ 



sin w p [n + + (1 - 4)S>] + cos w p [X* + X p + (4e 2 p + l)S eh ] , (A3) 



[Z* + Z p + Z GR + (1 - e 2 p )(AS ee - S qq )] - cosu; p —- [Y* + Y p + (l- e 2 p )S qh ] 

(XT Sill J. 

- sintup^r [X* + X P + (4e 2 + l)S eh ] , (A4) 
§ = ^[y* + Y P + (l-el)S qh ] + S ^[X^ + X p + (Ael + l)S eh ], (A5) 
^2 = ^-y^ + X^ + W^), (A6) 

= |(-^ + ^7 + W (A7) 

In eq. (A6) - (A7), the subscript 7 runs through the three spatial indexes, E, Q, and H, with E, Q, H being the three axis 
defining the orbital plane of the binary companion, and e, q and h defining that of the planet orbit, and = e • E, etc. 
The inclination / is taken to be that between the two planes, with the binary plane held constant during the evolution. 
The other two angles, oj p and f2, are the argument of the planet pericentre, and the longitude of the planet ascending 
node, respectively. The spin vector ft* (O p ) belongs to the host star (planet), and E = r g ±M±R 2 is the stellar moment 
of inertia, with r g + typically called the gyro-radius. I p is for the planet. The planet orbit has an angular momentum of 
magnitude h p = [G(M* + M p )a p (l — e^)] 1 / 2 , and the reduced mass of the planet-star system is fi = M*M p /(M* + M p ). 
All other symbols are defined in EKE. We do not repeat these definitions here. 

We introduce two commonly used dimensionless numbers, ki (the tidal Love number) and Q (the tidal quality factor), 
and relate them to the notation in EKE, 







A* 


= k 2 *Rl 
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kzpRp 
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3fc 2p 




Q* 


\a p/ 




t Fp 


Qp 



(A8) 



Rp\ 5 M* 



M„ »>■ (A9) 



DYNAMICAL TIDE 



We provide a simple description of the dynamical tide treatment. The final results are presented in §3.3. 

The period of the stellar spin is likely long compared to both the sound crossing time and the dominant tidal forcing 
period. So we assume a non-rotating star. Following Press & Teukolsky (1977) and Lai (1997) and notations therein, we 
decompose the tidal potential at a point r = [r, 6, </>] inside the star, due to a planet at location D(t) = [D(t), tt/2, f(t)], 
as 

£ 

U(r,t) = -GMpJ^Wem^YUd^)^- 1 " 1 ^ . (Bl) 
The stellar fluid responds with a displacement £(r, i) which can be projected onto various free eigen-oscillations, 

£(r,t)=£M*)W r )- (B2) 

a 

The eigenfunctions are normalized as J d 3 xp^ a ■ £* = 1. We denote the mode eigenfrequency as a a and damping rate as 
7 Q , and introduce a a (a) which is related to a a (t) by 



L 



a a (t)= / dercxp~ t<Tt a a (a). (B3) 



Substituting equation (B2) into the fluid equation of motion, we obtain 



a a (a) = J + ^ £ ^QatKW), (B4) 



where the periapse distance D p = a p {\ — e p ), Ki m (a) is the frequency spectrum for the Y^ m component of the tidal 
potential, and Q a i is the spatial overlap integral between the tidal forcing and the a eigenfunction, 



dt 



i+i 

exp[— imf(t) + iat], (B5) 



Qai = J d 3 x pC a ■ V{r l Y, m ) = J d 3 x p2 r e Y em . (B6) 
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Here, p' a is the Eulerian density perturbation for mode a. The continuity equation gives V • (p£ a ) — —p' a - 

We first focus on Ki m . For a strictly periodic orbit, K^ m (a) equals a sum of delta functions with non-zero values only 
at the integer multiples (j) of the orbital mean motion. 



K im (a) = y^Cj(a) S((T-jn p ). 
3=0 

To obtain Cj, we calculate the tidal spectrum integrated over one orbital period, 



(B7) 



Win 
27T 



r-2rr/n p 
/ dt 

Jo 



D(t) 



e+i 



exp[— zm/(i) + iat]. 



(B8) 



Unlike Ke m , this quantity is a smooth function of frequency For £ = \m\ = 2, the forcing maximum lies around 
a ~ (1 — e p )~ 3 / 2 \m\n p , i.e., a ~ 2n p for e p <l and a ~ 90n p for e p = 0.927. This is shown in Figure 3. When integration 
in K( m is carried out over N number of orbital cycles, the resulting forcing spectrum can be approximately described by 
the following log-normal distribution, 



3 = 1 



(B9) 



And when N — > oo, we retrieve Ki m with Cj(a) = \pKn v K\ m {o) 

Substituting the above expression for Ki m {a) into equations (B2)- (B4), we obtain the mode amplitude, 



C3+i)np 6{<j-jn p )e 
da — 



GM, 



(J- 



2ij a a Dl +1 



jV^npQ a eKj m (aj 



(BIO) 



where Uj = jn p is the closest orbital multiples to the mode frequency a a . 

The value of the tidal overlap, Q a £, decays sharply with both I and the radial order (n) of the mode. For a high order 
gravity-mode, most of the net contribution to Q a i comes from the upper evanescent region of the mode - the bottom of 
the surface convection zone where p' / p is fairly constant. We find that Q a f, scales with n as 



Qai oc n 



(BH) 



